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We investigate the entanglement between any two spins in a one dimensional Heisenberg chain as 
a function of temperature and the external magnetic field. We find that the entanglement in an anti- 
ferromagnetic chain can be increased by increasing the temperature or the external field. Increasing 
the field can also create entanglement between otherwise disentangled spins. This entanglement can 
be confirmed by testing Bell's inequalities involving any two spins in the solid. 
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It is well known that distinct quantum systems can be 
correlated in a "stronger than classical" manner [H-p|. 
This "excess correlation", called entanglement, has re- 
cently become an important resource in quantum infor- 
mation processing Q| . Like energy, it is quantifiable [g|-0 . 
This motivates us to ask how much entanglement exists 
in a realistic system such as a solid (the likely final arena 
for quantum computing |8|]) at a finite temperature. The 
ID Heisenberg model [pi W| is a simple but realistic |0] 
and extensively studied |p^-|l5[ solid state system. We 
analyze the dependence of entanglement in this system 
on temperature and external field. We find that the en- 
tanglement between two spins in an antiferromagnetic 
solid can be increased by increasing the temperature or 
the external field. Increasing the field to a certain value 
can also create entanglement between otherwise disen- 
tangled spins. We show that the presence entanglement 
can be confirmed by observing the violation of Bell's in- 
equalities. However, on exceeding a critical value of the 
field, the entanglement vanishes at zero temperature and 
decays off at a finite temperature. In the ferromagnetic 
solid, on the other hand, entanglement is always absent. 
We compare the entanglement in these systems to the 
total correlations. 

The entanglement of formation |q] is a computable en- 
tanglement measure for two spin-^ systems (qubits) [nw. 
We will use this measure to compute the entanglement 
between different spins in the ID isotropic spin-2 Heisen- 
berg model. This model describes a system of an arbi- 
trary number of linearly arranged spins, each interacting 
only with its nearest neighbors. Recently, entanglement 
in linear arrays of qubits have attracted interest 18-2(]|] 
and in Ref. [|9| the entanglement in the ground state of a 
Heisenberg antiferromagnct has been computed. But en- 
tanglement in the natural state of a system as a function 
of its temperature remains to be studied and the pos- 
sibilities of increasing this entanglement by an external 
magnetic field remains to be explored. The Hamiltonian 
for the ID Heisenberg chain in a constant external mag- 
netic field B, is given by 



where ct* = {al.,ay, a],) in which crl-iyu are the Pauli ma- 
trices for the ith spin (we assume cyclic boundary condi- 
tions l + iV=l). J>0 and J < correspond to the an- 
tiferromagnetic and the ferromagnetic cases respectively. 
The state of the above system at thermal equilibrium 
(temperature T) is p(T) = er"^!^^ jZ where Z is the par- 
tition function and k is Boltzmann's constant. To find 
the entanglement between any two qubits in the chain, 
the reduced density matrix p^ (T) of those two qubits is 
obtained by tracing out the state of the other qubits from 
p(r). Entanglement is then computed from p^{T) follow- 
ing Ref. ||l^. As p{T) is a thermal state, we refer to this 
kind of entanglement as thermal entanglement. Thermal 
entanglement is expected to behave differently from the 
usual solid state quantities (magnetization, correlations 
etc.), as the entanglement of a mixture of states is of- 
ten less than and at most equal to the average of the 
entanglement of these states. 




N 



H 



Y,{Bcjl + Ja\a 



■i+i\ 



(1) 



FIG. 1. We have plotted the entanglement E between two 
qubits interacting according to the antiferromagnetic Heisen- 
berg model as a function of the external field B and tem- 
perature (multiplied by the Boltzmann's constant) kT with 
coupling J = 1. The B — 4.6 line pointed out in the figure 
shows that for certain values of B it is possible to increase E 
by increasing T. AtT = 0,E has a sharp transition from I to 
as i? crosses the critical value oi Be — 4. E always becomes 
zero for values of T exceeding Tc = 8/k In 3 



We first examine the 2 qubit antiferromagnetic chain. 
We will use the entanglement of formation p, |l6|jl7| ] to 
calculate the entanglement of the two qubits. To calcu- 
late this entanglement measure, starting from the density 
matrix p, we first need to define the product matrix R of 
the density matrix and its time-reversed matrix 



R = p{(Jy ® ay)p*{<Jy ® U y) . 

Now concurrence is defined by 

C ~ max{\i — A2 — A3 — A4, 0} 



(2) 



(3) 



where the A^ are the square roots of the eigenvalues of R, 
in decreasing order. In this method the standard basis, 
{|00), |01), |10), |11)} must be used. The entanglement of 
formation is a strictly increasing function of concurrence, 
thus there is a one-to-one correspondence. The amount 
of entanglement in our special case is given by 
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where C is the concurrence given by 
C = if e^'"'^ < 3, 
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Fig.n shows the plot of this entanglement as a function 
of magnetic field and temperature. 





FIG. 2. The topmost plot shows the variation of nearest 
neighbor entanglement E with B for A'^ = 6, kT = 0.1 and 
J = 1. The middle and the bottom most plot show the same 
for next nearest and next to next nearest neighbors respec- 
tively. The reason for the shapes of the curves is presented in 
the text. Note that Ie is 1 lattice spacing for all values of B 
below Be =3.24 and changes to 3 lattice spacings for a range 
of B after Be ■ This means one can magnetically tune in the 
entanglement between any two qubits by increasing B. We 
also see the decay of all types of entanglement shortly after 
B = B, = 4. 




FIG. 3. We have plotted the variation of next to nearest 
neighbor entanglement E with B at kT = 0.1 and J — 1 for 
three values of A'^. We see that greater the A'^, a larger value 
of B is needed to tune in this entanglement, which is absent 
until B reaches a certain value. However, this entanglement 
disappears irrespective of A'' shortly after B exceeds Be = 4. 



For B = 0, the singlet is the ground state and the 
triplets are the degenerate excited states. In this case, 
the maximum entanglement is at T = and it decreases 
with T due to mixing of the triplets with the singlet. For 
a higher value of B, however, the triplet states split, and 
1 00) becomes the ground state. In that case there is no 
entanglement at T = 0, but increasing T increases en- 
tanglement by bringing in some singlet component into 
the mixture. On the other hand, as B is increased at 
T = 0, the entanglement vanishes suddenly as B crosses 
a critical value of Be = 4 J when |00) becomes the ground 
state. This special point T = 0, B = Be, at which en- 
tanglement undergoes a sudden change with variation of 
B, is the point of a quantum phase transition |^l[ (phase 
transitions taking place at zero temperature due to vari- 
ation of interaction terms in the Hamiltonian of a sys- 
tem). At any finite T, however, entanglement decays 
off analytically after B crosses Be- In the ferromagnetic 
case, the state of the system at _B = and T = is 
an equal mixture of the three triplet states. This state 
is disentangled m. Increasing B increases the propor- 
tion of 1 00) in the state which cannot make it entangled. 
Increasing T increases the proportion of singlet in the 
state which can only decrease entanglement by mixing 
with the triplet. Thus we never find any entanglement in 
the 2— qubit ferromagnet. These features of the 2— qubit 
Heisenberg model are also present in the N qubit model 
(which we investigate numerically) along with additional 
features, which we describe next. 

We first plot (Fig.g) how the entanglement between 
nearest, next nearest and next to next nearest neigh- 
bors in an antiferromagnet vary with B for a finite but 
low T (so that the entanglement is predominantly deter- 
mined by the ground state). For the nearest neighbor 



entanglement there are dips in the entanglement at cer- 
tain points. These dips are due to the mixing of two 
different entangled ground states at these points. After 
exceeding a certain value of B (say, -Be, which might 
depend on N), an equal superposition of states with 
only one spin up becomes the ground state. This state 
l^sym) = ^(|100...0) + |010...0) + ... + |000...1)) has 
entanglement between any two pairs. Thus we see the 
next nearest and the next to next nearest neighbor en- 
tanglement becoming finite only after B crosses Be- One 
can call this entanglement between non nearest neighbors 
magnetic entanglement as it is brought about by increas- 
ing B. When B is increased further, beyond a critical 
value Be = 2 J{1 + i±tl)^ + ^^^^^^ cos ir/N} < 4 J the 
disentangled state |00...0) becomes the ground state. At 
precisely T = 0, crossing Be ensures the complete van- 
ishing of all types of entanglement. For finite T, all types 
of entanglement decay to zero gradually after Be- This 
is illustrated in Fig.^ An interesting point, shown by all 
our numerical evidence, is that the change in entangle- 
ment AE at constant temperature due to change in B, 
can never exceed \AB\/kT. This might not be surprising 
as entanglement is an entropic quantity and |A_B| is the 
change in internal energy. 




FIG. 4. The nearest neighbor entanglement E is plotted as 
a function of kT at_B = 4.2, J = l and various values of A*' 
(from top to bottom, iV = 6, TV = 8, A^ = 10, Af = 9, iV = 7 
and A*' = 5). This graph shows that E of even A'^ states de- 
creases, E of odd A' states increases with A'^ and they both 
tend to merge with each other for high A'' (they almost coin- 
cide for A'^ = 9 and A'^ = 10). The plot also illustrates that 
one can increase the entanglement by increasing T. 



If we define a quantity called the entanglem,ent length 
Ie as the smallest separation between qubits beyond 
which the entanglement disappears, then for a small 
range of B after crossing Be, Ie becomes equal to the far- 
thest neighbor separation (i.e it can be made arbitrarily 
large). We have checked this numerically upto N = 13, 
and it is reasonable to conjecture that this will be true 



for any N . If this conjecture is false, it will still be inter- 
esting to find the value of N beyond which you can never 
increase Ie to the largest neighbor separation. Of course, 
as evident from Fig||, the further the qubits are, lesser is 
the magnitude of the entanglement between them. Note 
that the above definition of entanglement length differs 
from that defined in Ref. p2| where quantum to classical 
transitions in noisy quantum computers was studied (see 
also Ref. [p3[ for transitions in quantum networks). 




FIG. 5. This graph shows the variation of total mutual in- 
formation with temperature for the antiferromagnetic (AF,I) 
and ferromagnetic (F,I) case for A' = 10, B = 4.2 and \J\ — 1. 
The entanglement for the antiferromagnetic case (AF, E) is 
also plotted for a comparison. 



As mentioned earlier, l^sym) becomes the ground 
state for a certain range of values of the external mag- 
netic field (confirmed numerically up to A'^ = 13 and 
conjectured for other values of N). At extremely low 
temperatures and appropriate magnetic fields, thus, the 



state of the chain will almost be 1^; 



sym/ 



This state has 



the interesting property that there exists entanglement 
between any two qubits. The reduced density matrix of 
any two spins in the state |^sym) is p = ■^|^+)(^+| + 
(l-|r)|00)(00| where |*+) = ^(|01) + |10)). This is an 
entangled state. For any N, if we measure the state of all 
qubits except two, those two qubits would be projected 
onto a maximally entangled state, which can then be ver- 



24]. Even for any 
y or magnetically 



ified through Bell-CHSH inequalities | 
other mixed state which may thermal 
generated, there exists A neccessary and sufficient con- 
dition to check whether the CHSH inequality is violated 
P4| . Of course, one has to make an appropriate choice 
of measurement axes on the two spins in the solid. As 
different components of the magnetic susceptibility ten- 
sor are proportional to spin-spin correlations in different 
pairs of directions [^ , a CHSH inequality can tested by 
measuring different components of the magnetic suscep- 
tibility tensor . 

We now look at the dependence of entanglement on 
T in the N qubit case for a fixed B. Figjj shows that 



one can increase entanglement by increasing T. After a 
certain T, all entanglement dies out. In all simulations 
we find this temperature to be lower than Tc — 8/ In 3. 
Also, we see that the curves for entanglement in the case 
of even and odd N approach each other as N increases. 
This seems reasonable because for large N, it should not 
make a difference to the nearest neighbor entanglement 
whether we add or subtract a qubit somewhere far in the 
chain. As with the 2-qubit case, we find no thermal or 
magnetic entanglement in a ferromagnetic chain. 

We would now like to compare the amount of en- 
tanglement in the solid to the total two qubit correla- 
tions. An information theoretic measure of these cor- 
relations is the mutual information given by I{i : j) — 
S{p'^) + S{p') — S{p^^) where p*, p' are the density ma- 
trices of the ith and the jth spin respectively, p^^ is their 
joint state and S{p) represents the entropy of p. In 
a manner similar to the connected correlation function 
pa, this quantity measures the effect that genuinely re- 
sults from the interaction between particles. A plot of 
I{i : j) with temperature is shown in Fig.^. It is inter- 
esting to note that though entanglement is always ab- 
sent in a ferromagnet, / for nearest neighbors (stemming 
entirely from classical correlations) can be increased by 
increasing the temperature. It is well known that the 
magnetic susceptibility is proportional to the spin-spin 
correlations |l^]. It would be interesting to investigate 
whether any difference arises between the antiferromag- 
netic and ferromagnetic susceptibility tensors due to the 
complete absence of entanglement in the latter case. 

In this letter, we have introduced the concepts of ther- 
mal and magnetic entanglement and analyzed their be- 
haviour in the ID isotropic Heisenberg model. We have 
found critical values of field beyond which entanglement 
disappears at zero temperature and declines at finite tem- 
perature. Our results indicate that there is also a crit- 
ical temperature after which all entanglement vanishes, 
though there is a range of field in which entanglement can 
be increased by increasing the temperature. Based on nu- 
merical evidence, we have conjectured that the entangle- 
ment length can be made arbitrarily large by applying an 
appropriate external magnetic field. We have also com- 
pared the total correlations to the entanglement. Our 
work raises a number of interesting questions and conjec- 
tures to prove and the possibility of numerous general- 
izations such as higher dimensions, non nearest neighbor 
interactions, anisotropics, other Hamiltonians and so on. 
In addition, we also showed that by applying a suitable 
magnetic field and lowering the temperature sufRciently, 
and doing suitable projections, one can create a state 
which violates the Bell's inequalities. The "natural" en- 
tanglement can be verified in these cases. In future we 
will investigate how to map this natural entanglement 
out of the solid (eg. by neutron scattering) and use it as 
a resource in communications. 
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